An Area Law for One Dimensional Quantum Systems 



M. B. Hastings 
Center for Nonlinear Studies and Theoretical Division, 
Los Alamos National Laboratory, Los Alamos, NM, 87545 

We prove an area law for the entanglement entropy in gapped one dimensional quantum systems. 
The bound on the entropy grows surprisingly rapidly with the correlation length; we discuss this 
in terms of properties of quantum expanders and present a conjecture on completely positive maps 
which may provide an alternate way of arriving at an area law. We also show that, for gapped, 
local systems, the bound on Von Neumann entropy implies a bound on Renyi entropy for sufficiently 
large a < 1 and implies the ability to approximate the ground state by a matrix product state. 

66 '■ 

There are many reasons to believe that the entanglement entropy of a quantum system with a gap obeys an area 
law: that the entanglement entropy of a given region scales as the boundary area, rather than as the volume. In one 
dimension, conformal field theory calculations show that away from the critical point the entanglement entropy is 
bounded, diverging proportionally to the correlation length as a critical point is approached 1]. In higher dimensions, 
systems represented by matrix product states 0,0] obey an area law. 

However, despite this, there is no general proof of an area law. This is somewhat surprising, since it has been 
, proven that correlation functions in a gapped system decay exponentially^, and one might guess that the decay of 
correlation functions implies that only degrees of freedom near the boundary of the region may entangle with those 
outside. However, the existence of data hiding states Q shows that one can have states on bipartite systems with 
O ,■ small correlations and large entanglement. Further, the existence of quantum expanders^, 0] shows that one may 
have matrix product states in one dimension that have all correlation functions decaying exponentially in the distance 
between the operators, but still have large entanglement. This indicates some of the difficulty in proving an area law. 

At the same time, for a gapped system to violate an area law would require some very strange properties. For one 
thing|3], the thermal density matrix can be well approximated by a matrix product operator. This implies that unless 
a plausible assumption 6| on the density of low energy states is violated, the ground state can be well approximated 
by a matrix product state. 

In this paper, we succeed in providing a proof of an area law for one dimensional systems under the assumption of a 

■ gap. The result, however, bounds the entanglement entropy by a quantity that grows exponentially in the correlation 
p\| ", length. This is much faster than the linear growth one might have expected. We will comment later on why this 

bound might in fact be reasonably tight. In the process of deriving this result, we will derive bounds for gapped local 
systems which inter-relate three quantities: the Von Neumann entropy, the Renyi entropy, and the error involved in 

■ approximating the ground state by a matrix product state. 

We begin by defining the lattice and Hamiltonian. We consider finite range Hamiltonians for simplicity. It is likely 
that the results can be extended to exponentially decaying interactions, but for simplicity we do not consider this 
here. In fact, having decided to consider only finite range interactions, we may group several sites into a single site, 
and thus simplify to a problem with only nearest neighbor interactions. 

Specifically, we consider a finite volume one-dimensional lattice, with sites labeled i = 1,2, ...,N, with a D- 
dimensional Hilbert space on each site. We consider finite-range Hamiltonians of the form H = The 
finite range condition is that -f/i^+i has support on the set of sites i and i + 1. We additionally impose a finite 
interaction strength condition bounding the operator norm, that the operator norm H-ff^+il < J for some J. 

The properties imply a Lieb- Robinson boundjijl. [Tot 11 1: there exists a velocity v and length scale £c such that for 
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any two operators A, B with support on sets X, Y respectively, 

|p(t),B]|| <cx |X|p||||B||expK c dist(X,Y)] (1) 

for \t\ < l/v where the distance between sets X,Y = min^x, jeY {\i — j\), where c is a numeric constant of order unity, 
and where A(t) = exp[iHt]Aexp[—iHt\. The velocity v will be of order J, while £c will be of order unity. 

We now introduce some notation. We let Xj± denote the set of sites i with j < i < k. We let denote the ground 
state of the Hamiltonian H and we let p® N = Wo) be the ground state density matrix. We let p° k denote the 
reduced ground state density matrix on the interval Xj^. That is, p°- k = tr^Xj k (Pi jv)> where the partial trace is 
over sites not in Xj^. We define the entropy of any density matrix pj k by S(pj,k) — t^ieXj fc (/°j,fc in (/°j,fe))- 

Theorem. Consider a Hamiltonian satisfying the finite range and finite interaction strength conditions above. Sup- 
pose H has a unigue ground state with a gap AE to the first excited state. Then, for any i, 

Sipli) < S max (2) 
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where we define 

S max = c eHi')MD)2t'^ D \ (3) 
for some numerical constant Cq of order unity, and where we define 

£ = max(2i;/A£,fc), (4) 

e = 

We do not consider the case of a degenerate ground state, but we expect that the theorem can be strengthened to 
include this case also. 

PROOF OF MAIN THEOREM 

To prove the theorem, we assume that it is false, so for some we have S(pi f ) > S rnax . Then, for any k > i we 
have S(p® k ) > S max — (k — i) ln(_D), and therefore, for all k with i < k < i + Iq where 

lo = S max /3\n(D), (5) 

we have 

S(p° 1<k ) > 2S max /3 = S cut . (6) 

Then, define Si, for each / < Iq, to be the maximum entropy of an interval of length I contained in Xi^ + i . That is, 
Si is the maximum of S(p® +1 over j such that Xj + i^ + i C Xij + i Q . Clearly, 

Si < HD). (7) 

Further, for any I, we have 

S 2 i < 2S t . (8) 

However, it should be apparent that it is not possible for Eq. © to be saturated. If S21 = Si, then we have for some 
j that S(j)j_i +1 = S(j)j_i +1 j) + S(pj +1 j + i)- In this case, the density matrix P]~i+ij+i is equal to the product 
of density matrices /0°_ i+1 j ® p® +1 This implies that tv{Hp\^ £g> p°- +1 N ) — tv(Hp\ N ) since H is a sum of terms 
Hi t i + i acting on nearest neighbors. However, given a unique ground state, this implies that n = Pi,j ® Pj+i jv> 
contradicting the assumption of a non- vanishing entanglement entropy S(pij). The key idea of the proof will be to 
improve on Eq. ((5|) even further, and use the fact that S(p1 •) is not only non- vanishing, but bounded below by S cu t 
to show that the following stronger claim holds for all I <Iq: 

S 2 i<2S l -(l-2Ci(0exp(-l/e))l/e + HCi(0) + C2, (9) 

for some function Ci(£) which is bounded by a polynomial in £, and some numerical constant C 2 of order unity. 

Once Eq. (J9j) is shown, we can iterate it. The dominant terms on the right-hand side of Eq. ([9]) at large I are 2Si~l/£', 
so that one can guess from the behavior of these dominant terms that we will find Si w I ln(D) — log 2 (l / £o)l / £' , up to 
subleading terms. Indeed this is the case. Set £0 = £'2 ln(2Ci(£)). We use Si < l\n(D) for Z = £0 as an initial bound 
and iterate to I = 2£ , 1 = 4£ , and so on. Note that 2Ci(£) exp[-£ /£' < 1, and Y^Lo 2C i(0 ex P[-2"W?']£o/£' 
is bounded by 1 + e~ 2 + e~ 4 + ... < 2. It follows that for I > £ , 

Si<\n(D)l-l[\og 2 (l/^)\/^ + (2 + ln(C 1 (0)+C 2 )l/^- (10) 

However, since Si must be positive, a contradiction will arise when U°g2G/£o)J > ^ n {D)£,' + (2 + C 2 )^' + 1/2 and 
so 

l < e In(2Ci(0)2^' m(o)+(2+c 2 )e'/fo+i/21 _ (u) 

Combined with Eq. ([5|) this gives the main result for some constant cq. 

Thus, we must now show Eq. To show this, it suffices to show that for all j, I with Xj+\ j+i C Xij + i that, 

S( P U+i, j+ l) < S(pU+i,j) + S(P° j+ i, j+ i) - (1 - 2C X (0 cM-l/m/e + HCi(0) + C 2 . (12) 
At this point, we need two lemmas that rely on the Lieb-Robinson bound and the existence of a gap. 
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Lemma 1. Let H be a Hamiltonian that satisfies the finite range, finite interaction strength, and gap conditions 
above. Then for any j and any I there exist Hermitian, positive definite, operators Ob(j, Z), Ol(j, Z), (9r(j, /) with the 
following properties. First, \\0 B (j,l)\\ < 1, \\0 L (j,l)\\ < I, and \\0 R (j,l)\\ < 1. 
Second, 

\\O B (j,l)O L (j,l)O R (j,l)-P \\ < dCOexpI-Z/e'] = e(0, (13) 

where Pq = ^oK^o * s the projection operator onto the ground state and the function Ci(£) is bounded by a polynomial 
in£. Finally, Ol{],1) is supported on X± t j, On(j,l) is supported on Xj + \ jq , andOB{j,l) is supported on Xj_; + i J+ ;. 
We prove this lemma in the appendix. 

The next lemma provides a way of bounding the entanglement entropy of a given region, given a certain assumption: 
that one can approximate the ground state to a certain accuracy in Hilbert space norm by a state with given Schmidt 
rank. In general, if one knows that a matrix product state provides a good approximation to the ground state, this 
lemma can be used to "bootstrap" that result into a bound on the entanglement entropy of the ground state. 

Lemma 2. Let 

k k 

P = J2 P (T) 7 ) <g> V R (a, MP, 7)*t(A l) ® ^O 9 . 7) (14) 

7 a=l 0=1 

be some density matrix with unit trace. Here, ^i(a,7) are states on X\.j and ^ R (a,y) are states on Xj + i,^. Then, 
we say that p is a mixture of pure states with Schmidt rank at most k. Suppose that {9q,p^!q) = P > 0, where \l/o *s 
the ground state of a Hamiltonian that satisfies the finite range, finite interaction length, and gap conditions above. 
Then, 

S(p° hj ) < \n(k)+Z>\n(2C 1 (0 2 /P)HD)+F(t;',D), (15) 

where 

F(£, D) = (£' + 4) \n(D) + 1 + ln{D 2 - 1) + m(£'/2 + 1), (16) 

and a similar bound holds for Renyi entropies with sufficiently large a as discussed below (here, a refers to the order 
of the Renyi entropy, not to a particular state; the particular use of a should be clear in context). 

Proof. To prove this, note that for any m the positive definite operator p{rn) = 
ObU, "m)Oh{j, ti)Oh(j, m)pO R (j, m)Oi(j, m)Os(j, m) has the following properties. First, p{m) is a mixture 
of pure states with Schmidt rank at most kD 2m . Second, from Eq. (flB"]) , tr(P p( m )) > (V? — e (0) 2 an d 
tr((l - P )p(m)) < e 2 and thus 

(* ,p(m),* ) ^ 1 ^ _ 2 



tr(p(m)) 



> l-2C 1 (e) 2 exp[-2m/^']/^ (17) 



Let the ground state \I/o be equal to J2 a =i Ao( a )^ L fi( a ) ® ^/R,o(ce), where the states ^L,o(ct) and ^>R,o(a) are 
again states on X\j and Xj + \^ respectively, with (^L,o(ct), *l,o(/3)) = (*fl,o(a), *H,o(/3)) = 5 a ,p- We order the 
different states such that if a < (3 then |A (a)| > |A (/3)| and we normalize so that J2 a \M{ a )\ 2 = 1- It follows from 
Eq. (JXTJ) that for all integer m 

\Ao(a)\ 2 <2C 1 (0 2 eM-2m/e]/P (18) 

a>kD 2m + l 



Let m! be the smallest integer m such that 2exp[— 2m/£']Ci(£) 2 /-P < 1. Thus, for m > 



in 



\Ma)\ 2 <exp[-2(m-m')/£']- (19) 

u>kD 2m + l 

We now maximize the entropy S(p® +; ) = — J2 a =i \Ao{ a )\ 2 l n (l^4-o( a )| 2 ) subject to the constraint (TIT))) . The 



maximum occurs when Yla=i l^o( a )| 2 = (1 — ex p[^2/£']) and for m > m! ', X)a=fcD 2m +i |^4o(ck)| 2 = (1 
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exp[— 2/£']) exp[— 2(m — m' )/£'], giving an entropy bounded by 

oo 

ln(fc) + (2m' + 2) ln(L>) ^(2(n - 1) ln(D) + ln(P> 2 - 1) + 2n/g) exp[-2n/^ / ](l - exp[-2/£']) (20) 

n=l 

-ln(l-cxp[-2/£']) 

= ln(fc) + (2m' + 2 + J^|ti0L) ln(£>) + + HD 2 - 1) exp[-2/£'] - ln(l - exp[-2/£']) 

< ln(fc) + (£' ln(2C x (0 2 /P) + £' + 4) ln(£>) + 1 + ln(D 2 - 1) + ln(£'/2 + 1), 

where we used the inequalities 2 exp[-4/£']/(l - exp[-2/£']) < (2/f ) exp[-2/^']/( 1 - cxp[-2/£']) < 1, and 
1/(1 - exp[-2/f']) < i'/2 + 1 for £' > giving Eq. (JTSJ) as claimed. 

We note that this proof can be extended to Renyi entropies S a (p\j)) defined by S a (px,j) = (1 — a) -1 ln(tr(p" ■)) 
for sufficiently large a. Maximizing the Renyi entropy subject to the constraint (fl~8|) gives S a {p\ ,•) < (1 — 
a)- 1 In{£^ (fcZ) w + 2 .D 2n ) 1 - O! (l - exp[-2/C']) Q exp[-2n/C'] Q }. The sum converges so long as 

cxp[-2a/£'](.D 2 ) (1 - Q) < 1, (21) 



in which case we have a bound on the Renyi entropy which differs from the bound bound (|15[) on the von Neumann 
entropy only in that the function P(£', D) is replaced by an a-dependent function. □ 

We now return to proving the main theorem. From Eq. (fT3|) it follows that (^o, Ob{j, Z)Ol(j, ^)Or(j, O^o) > l~ e (0 
and hence (* 0> B (j, Z)* > x (y O L (j,l)O R (j,l)V ) > (1 - e(Z)) 2 . Therefore, (¥„, B (j, l)*o) > 1 - 2e(i) and 
(*o,OlO'.OOk(J.O*o) >l-2e(i). Thus, 

tr(p°_ i+lj+i O B (j,Z))>l-2e(Z). (22) 

Set P = tr(P p° ■ <g> p° j+1 N ). From Eqs. (IHTTK)) . we find that £' ln(2Ci(£) 2 /P) ln(D) + P(£', D) > S cut . Therefore, 
P < 2d(0 2 expH5 cnf - F(e,D))/ln(D)C]. Let 

x = tr{0 B (j,l)pl j ® p° j+hN ) (23) 

and y = tt(p L (j,l)Oa (jM d ® P?+i,iv) > 1 - 2e(i). Then, P > tr(0 B (j, 0^0'. OOjtCfc OPi.j ® " e (0 > 

— \J x — x 1 \J y ~ y 1 — e(Z) as follows from a Cauchy-Schwarz inequality for the expectation value of a product of 
operators. Thus, P > x(l - 2e(Z)) - \/x - x 2 y/2e(l) - e(Z) > x(l - 2e(l)) - y/xy/2e{l) - e(l). Thus, 



x < {2C 1 (0 2 exp[-(S c „ t - P(£', £>))/ ln(£>)£'] + vWMO + 2e(Z)}/(l - 2e(Z). (24) 

Thus, Eqs. (|22|23p imply that the operator Os(j,l) has a large expectation value for the state . +J and a 

small expectation value for the state pj_ l+ltj <8> Then, the Lindblad-Uhlmann theorem[12j provides a lower 

bound on the relative entropy S(pV_ l+l j+l \\p?_ l+1J <g> p° +lj+; ) giving 

) (25) 
> (1 - 2e(Z)) ln((l - 2e{l))/x) + 2e(Z) ln(2e(Z)/(l - x)). 

Eq. (|25|) is the key step. Everything that follows consists of picking the constant C2 correctly. We first assume 
that lo < {Scut — P(£', D))/ ln(£>) — £'ln(Ci(£)). If this assumption fails, then since lo = S cu t/2\n(D), we have 
S cut /2\n(D) < F(€',D)/ln(D) + £'ln(Ci(£)), and thus by picking the constant c large enough Eq. (J3J) will still 
hold. Then, since Z < Z it foll ows fro m this assumption that 2Ci(£) 2 exp[-(S cnt - F(£' ,D))/ ln(D)£'] < 2e(/). Then 
from Eq. (24]), a; < (4e(Z) + ^/2e(Z)x)/(l - 2e(x)). Thus, as Z becomes large, we find that x and e(Z) approach zero 
exponentially. Thus, for large enough Z, we have S(p ( j_ l+1 j) + S(p® +1 j +l ) — S(p ( j_ l+1 j +l ) > (1 — 2e(Z)) ln(l/x) plus 
some constant of order unity and hence for all I we have 

S(pU+i,j) + S (P°+U+l) - S(P°-i+i^+i) > (1 - 2e(Z))ln(l/e(Z)) - C 2 , (26) 
where C2 is a numeric constant of order unity. This shows Eq. (|12[) and completes the proof. 
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MATRIX PRODUCT STATES 

Our main result is the bound ([2]). We now use the existence of a bound on entropy to construct an approxima- 
tion to the ground state by a matrix product state. At first, this might seem difficult, given that previous such 
constructions [HJ relied on the existence of a bound on the Renyi entropy, and we have a bound on the von Neumann 
entropy. However ideas similar to those used in the bootstrap lemma ([2]) will let us avoid these difficulties. 

The construction applies to any system given the finite interaction range, finite interaction strength, and gap 
conditions, and given the existence of a bound S(p1 j) < S max for all j. Consider a given j and write the ground state 
as in lemma ([2j by *o = J2a=i ^( a )^ L,o( a ) ® ^R,o{c), with *i,o(Q),tH,o(") being states on Xij and X j+ i^ N 
respectively and the states ordered so that |A (a)| is a non-increasing function of a. 

Suppose for some k' we have Ea=fe'+i \ A o( a )? > V 2 - The n, \A (k' + 1)| 2 < l/2k. Thus, S(p^) > (1/2) ln(2fc')- 
So, k' > exp(2S , ( j o° i ))/2. Hence, for 

fco =exp(25(p? J ))/2 (27) 

we have £ ° =1 |Ao(a)| 2 > 1/2. At this point we could directly apply lemma @ to bound the Renyi entropies, and 
thus get a matrix product form following results in[l3|], but we prefer to proceed more directly. We apply Eq. (fl8|) 
with P = 1/2, getting 

E |A («)| 2 <4C 1 (0 2 exph2m/r]. (28) 

a>k D 2m + l 

Therefore, 

£ \A (a)\ 2 < 4C 1 (0 2 exp(~Llog D (fc'Ao)J/e') <4C 1 (0 2 exp(l/0(fc , /fco) 1/C ' ln(D) . (29) 

a>k' 

This provides an estimate on how rapidly the Schmidt coefficients decay and hence how accurately the ground state 
may be approximated by a matrix product state. In particular, for a chain of length N, the error in approximating 
the ground state by a matrix product state of bond dimension k' scales as N(k' /fco) 1 ^ ln ^ D \ 

DISCUSSION AND A CONJECTURE ON COMPLETELY POSITIVE MAPS 

We now further explore the relationship between a gap, exponentially decaying correlations, and an area law. 
Consider a one dimensional system with a gap. We know that this state must have exponentially decaying correlations. 
However, this in itself does not imply an area law. For example, there exist matrix product states 

*(si,s 2 , — s N ) = £ A at 0(s 1 )A /3tl (s 2 )A 7i s(s 3 )... (30) 

a,/3,... 

with Si — 1...D in which the associated completely positive map forms a quantum expander so that a system with 
a low Hilbert space dimension D on each site and a short correlation length of order unity may have a very large 
entropy 0. Consider, however, the following model of a quantum expander modeled on that in [y]: we have a graph 
with coordination number D (which is a classical expander graph), where each node of the graph corresponds to a 
value of the bond variable a, each link from one node to another gets labeled with one particular value of s, and 
A a p(s) is non- vanishing only if the link from a to /3 is labeled with the given s. The total number of nodes in the 
graph is equal to the range of the bond variables, and we will denote this number by k. For a given value of a, 
the bond variable /3 can assume D different values, 7 can assume D(D — 1) different values, and so on, so that any 
correlations between a bond variable a and another far away bond variable which connects a distant pair of sites 
become small as required. This behavior is shown in Fig. 1 for a system with D = 3 and an arbitrary possible choice 
of a,/3,j,S,e. 

However, as we have seen in this paper, if there is entropy across a bond variable in a gapped, local system, then 
there must be mutual information of order £/£' between the set of sites within some distance I to the left of the bond 
and the set of sites within some distance I to the right of the bond. However, in the given expander map, there is 
no such mutual information. Consider a given set of sites i = i m iniimin + E ■■■,imax, with i max — i m in << N and 
_D* max ~*""" much smaller than the range k of the bond variable a. Trace out all sites outside to construct a reduced 
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FIG. 1: Illustration of an expander graph, with a particular choice of a,f3, ... marked. Note that each Greek index marks a site 
which is one site away from the previous index. 

density matrix on the given sites: the result is proportional to the identity matrix and any of the D lmax ~ lmin+1 states 
is possible with equal probability. Thus, this expander map is not the ground state of a Hamiltonian with an energy 
gap of order unity and an interaction range of order unity, even though it is the ground state of a gapped Hamiltonian 
with an interaction range of order \og D (k). In a sense, this expander map is mixing too quickly to be the ground state 
of a gapped Hamiltonian. 

There is, however, an interesting alternative way of seeing that this expander map cannot be the ground state of 
a gapped local Hamiltonian, which will lead to a conjecture we have on completely positive maps. We begin with a 
result on the decay of a certain kind of correlation function. 

Lemma 3. Let H a Hamiltonian H which satisfies the finite range, finite interaction strength, and gap conditions. 
Let be the ground state of H. Suppose = S a =i Ao(a)H! l,o{ol) <E> ^R t o{a), where ^l,o(o:) are orthonormal states 
on Xij and ^R,o( a ) are orthonormal states on Xj +1 ,N. Let B-l = J2 a =i O(a)^ l,o(®))(^ l,o ® Ar, where 1r is the 
unit operator on Xj + i,n, for some function 0(a). Similarly, let Br = 1l ® J2 a =i O(a) , I'fl,o(a0)(^ , -R,o- Suppose 
\\Bl\\ < 1; so that \0(a)\ < 1 for all a. Finally, let A be an operator with support on -X'lj-i U Xj + i + \^ and with 
\\A\\ < 1. Then 

(*o,A5 L f ) - (*o,^*o)(*o,B L *o> <3y/te(t) + e(l), (31) 
where e(l) is given as before by Ci(£) exp[— Z/£']. 

Proof. To prove this, dchne OB(J,l),Oi,(j,l),0 R (j,l) as in Lemma 1. Recall that (\I>o> Os{j, O^o) > 1 — 2e(Z) and 
(*o, R (j, 0*o) > 1 - 2e(0- Hence, \{0 B (j, I) - l)*o| < V 2e ( l ) and \(Or(j> - !)*o| < VM 1 )- We now use a series 
of triangle inequalities. First, 

| (* , AB L *o) - (*o, AB l O r (j, /)f o) I < VHl)- (32) 

Next, 

\(*o,AB L R (j,l)* ) - (*o,0 B (j,l)AB L R (j,l)* )\ < VW)- (33) 
However, [0 R {j,l),B L ] = and [0 B (j,l),A] = 0. Also, B L V = B R ^ . Thus, 

(*o, B {j, l)AB L R (j, 0*o> = (*o, AO B (j, 1)OrU, 1)Br*o)- (34) 

Next, 

\(*o,AObV,1)OrV,1)B r *o) - ^ ,AOb(3,1)Or(3,1)BrO l {3,1)* )\ < y/to®. (35) 
Using [B R , L (j, I)} = we find that 

\(*o,AB L * ) - (* ,AO B (j,l)O R (j,l)O L (j,l)B R * )\ < 3^(1). (36) 
Thus, \(^ ,AB L ^ ) - (AP B L ty ) \ < 3 A /2e(/) + e(l), completing the proof. □ 
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Now, this lemma also implies that the matrix product state arising from the given expander map cannot be the 
ground state of a Hamiltonian with gap and Lieb-Robinson velocity of order unity. In a slight abuse of notation, let 
a be the bond variable which joins site j — I to j — I + 1, let 7 be the bond variable which joins site j to j + 1, and 
let e be the bond variable which joins site j + I to j + I + 1. This is shown in Fig. 1 for a system with I = 2. 

For any given value of the bond variable a, the bond variable 7 can have a wide range of possible values, of order 
(D — 1)'. However, for given values of both a and e, the possible range of values of 7 is much more restricted, as 7 
will tend to lie on the shortest path joining a to e as shown in the figure. In this manner, it is possible to construct 
operators A, Bl such that Eq. (|3Tj) is violated for this state, although we omit the detailed construction. 

In fact, we have not been able to find any matrix product states which obey Eq. (j3Tj) while still having a large 
entanglement entropy. We thus make the following conjecture on completely positive maps: 

Conjecture 1. There exists a function f(D e ff) with the following property: consider a matrix product state ^0 o,s 
in Eq. H30\). Suppose that this state satisfies Eq. \31}) for all A, Bl for some given Then, the entropy S(pij) 

is less than or equal to f(D^ ). 

If this conjecture were shown, then it would give a different way to prove an area law. 



CONCLUSION 



In conclusion we have given a proof of an area law, although the bound is quite weak. We note that the bootstrap 
lemma ([2]) gives some intuitive idea as to why it is difficult to prove bounds on the entanglement entropy. Given a 
good approximation to the ground state with a state of given Schmidt rank, we have a bound on the ground state 
entropy. Also, given a bound on the ground state entropy, we can estimate how well we can approximate the ground 
state with a state of given Schmidt rank as in Eqs. (|27I29[) . However, this kind of argument leads to circular reasoning 
and thus does not help provide an area law. These arguments, do however, inter-relate the Von Neumann entropy, 
the Renyi entropy, and the error in approximating the ground state by a matrix product state in gapped system, as 
in Eqs. (|15|27)29|) . 

We finally consider the implications for numerical simulation of one dimensional quantum systems. If such a system 
has a gap, then the ground state is close to a matrix product state. However, finding the best matrix product state 
may be hard problem L|. Imagine, though, that we consider a family of Hamiltonians which start from a Hamiltonian 
with a known ground state and keep the gap open. Then, by following a quasi-adiabatic evolution [l5| along this path 
and truncating the matrix product state to keep the bond dimension bounded, we can well approximate the ground 
state at the end of the evolution 
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APPENDIX: APPROXIMATING THE PROJECTION OPERATOR 



In this appendix we show how to approximate the projection operator by a product Ob{j, 1)Ol(J, 1)Or(J, I) as 
described in Lemma 1. For simplicity, let us add a constant to H so that the ground state energy is equal to zero. We 
start by recalling the result in [15j . that it is possible, given a gapped Hamiltonian, to write the Hamiltonian as a sum 
of terms such each term approximately annihilates the ground state. We perform the derivation slightly differently to 
get a tighter bound (|39I40|) which depends only on the boundary of the terms. Let 

Hl = ^ Sm+Ij (37) 

i,i<j-l/3 

Hb — ^ Hij + i, 

i,j-l/3+l<i<j+l/3 
ia>j+l+l/3 
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so that H = Hl + Hb + Hr. We then choose to add constants to Hl, Hb, Hr so that (Hl) = {Hb) = (Hr) 
Define 



AE 

V^M J-c 
AE 



Hr 



fhrq 
AE 



dtH L (t) exp[-(tAE) 2 /2q], 
dtH B (t) exp[-(tAE) 2 /2q], 

1 
) 

dtH R (t) exp[-(tAE) 2 /2q], 



(38) 



where we pick q = (l/3)AE/(2v). Note that H°+H B + H R = H. Let b L = [H,H L \. Note that ||6 L || < J 2 . Then 

\Hl* \ < AE- 1 \HH° L y \ = AE- 1 \[H,Hl}y \ (39) 
AE 



= AE' 



f2jxq 



dtb L (t)^ \ 



< A^- 1 J 2 C(exp[-V3C]), 

where on the last inequality we used the assumption of a gap and 0(...) is used to denote a bound up to a numeric 
constant of order unity. Similarly, 



\H B ^ \ < 0(AE~ 1 J 2 exp[~l/3^}), 
\H B V \ < 0(AE- 1 J 2 exp[-l/3^}). 



(40) 



Using the Lieb- Robinson bound, for the given value of q it is possible to approximate , H B , H B by operators 
Ml, Mb, Mr respectively such that Ml — Hl is supported on Xj_ 2 i/z.j, Mb — Hb is supported on Xj_2//3,j+i+2Z/3! 
and M R - Hr is supported on X j+1 3+1+2l/3 and such that \\M L - H L \\ < 0(AE~ 1 J 2 exp[— Z/3£]), \\M B - H B \\ < 
OiAE-^Fexvl-l/Zt}), \\M R -H R \\<0(AE- l J 2 exp[-l/3^]), 

Thus, we have |M L * | < C(AE^ 1 J 2 exp[-Z/3£]), and similarly for M B and Mr. We now define L (j, Z) to project 
onto eigenvectors of Ml with eigenvalue less than or equal to AE~ X J 2 exp[— Z/6£], and define R (j, I) to project onto 
eigenvectors of Mr with eigenvalue less than or equal to AE -1 J 2 exp[— Z/6£]. Thus, 



\(0 L (j,l) -1)* \ < £>(exp[-Z/6£]), 
\(0 R (j,l) - l)*o| < 0(exp[-l/6£]). 

Clearly, Ol(j,1) and R {j,l) are supported as required by Lemma 1. 
We now define an approximation to the projection operator 



(41) 



P, 



AE 
■fUtq 



J dtcxp[i Ht] cxp[-(tAE) 2 /2q]. 



Using the spectral gap we have \\P q — Foil < 0(cxp[—l/3(,]). Thus, using \\Ml 
C(AF- 1 J 2 exp[-Z/3C]), we have 



Me 



M fi 



AE 



sPntq 
AE 



J dt exp[i(M L + M B + M R )t] exp[~(tAE) 2 /2q] - P a 



(42) 
H\\ < 

(43) 



dtexp^ J dt' exp(i(M L + M R )t')M B exp(-i(M L + M R )t') 

exp[i(M L + M R )t] exp[-(tAE) 2 /2q] - P 
< O {AE- 2 J 2 ^exp[~l 
where the exponential of the integral over t' is i'-ordered. Thus from Eq. (|4T 



AE 



rHxq 



dt exp 



i / dt' exp(i(M L + M R )t')M B exp(-i(M L + M R )t') 
L Jo 

exp[i(M L + M R )t] cxp[-{tAE) 2 /2q]0 L {j, 1)O r {j, I) - P a 
< 0(AE- 2 J 2 ^/q~exp[-l/3£] + exp[-Z/6£]). 



(44) 
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However, \\exp[i(M L +M R )t]0 L (j,l)0 R (j,l)-0 L (J,l)0 R (j,l)\\ < 2|t|AiT 1 J 2 0(exp[-Z/6£]). Combining this with 
Eq. {44]) we find that 



AE 



dt exp 



/2vrg 

cM-(t&E) 2 /2q]0 L (j, l)0 R (j, I) - Po 
< 0(AE- 2 J 2 y/qex.p[-l/6£]). 



i / dt' exp(i(M L + M R )t')M B exp(-i(M L + M R )t') 



(45) 



Consider the operator 
AE 



Pb 



J dtexp[i J dt'exp(i(M L +M R )t')M B exp(-i(M L + M R )t')]e-xp[-{tAE) 2 /2q}. 



(46) 



Using a Lieb-Robinson bound for Ml,M r (and noting that the difference Ml + M R — Hl — H R has support on 
Xj-2i/3,j+i+2l/3)-> we can approximate P R by an operator Ob{j,1) with support on Xj-i+x^+i such that \\P R — 
B (j, Oil < 0(AE-\J^qeM-l/^])- Thus, 



I|O s (j,0Ol(j, Z)Ofl(j, - -Poll < 0(AE- 2 J 2 ^lAE/2veM-l/H] = 0(y/l/i? exp[-l/6{]). (47) 
This completes the result. 
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